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Back Flow From Jet Flumes in Vacuum 

By Norman T. Grier 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 

ABSTRACT 

Approximate analytical calculations are performed to 
determine the particle density at surrounding points for- 
ward of a jet nozzle in vacuum due to the gas exhausted in 
the jet plume. The calculations are made for helium, argon 
and nitrogen. The back flow is determined from the number 
of particles that leave the continuum regime and cross the 
particular point in question without experiencing any in- 
terparticle collisions. 


INTRODUCTION 

For spacecraft accelerating or maneuvering in space, 
the back flow from the jet exhaust may be a large contri- 
butor to the ambient environment of the spaceship. In 
order to get an estimate of the density of the back flow, 
an approximate calculation is performed for points located 
in the plane of the nozzle exit. 


The gas in the plume is considered to be continuum and 
Maxwellian until the density decreases to the point where 
the mean free path is comparable to the distance that the 
particular mass of gas has travelled from the nozzle exit. 
Afterwards, the flow is assumed free molecular with no in- 
terparticle collisions occur ing . The back flow is deter- 
mined from the number of particles that leave the defined 
continuum regime and cross the particular point P(x,y, z) 
in question. This model is depicted in the sketch on the 
following page. 
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Since most of the Plume constant Knudsen surface (de- 
fined by Eqs. (5) and (7) below) is usually far from the 
jet exit plane, it is assumed that the flow has nearly 
reached its adiabatic limit and has a radial direction; 
that is, the streamlines diverge from a common source 
point located at the nozzle exit. The validity of these 
assumptions is discussed in Refs. 1, 2, and 3. The 
gases used are helium, nitrogen, and argon. In all the 
calculations the jet exit Mach number is 5. 

ANALYSIS 

In Refs. 1, 2, and 3 it is shown that the density 
along the plume axis of a jet exhausting into a vacuum 
is approximately given by 



where r* is the nozzle throat radius, r the distance 
from the source point, p c is the chamber density, and 
(Ref. 2) 



where 
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0^ is the thrust coefficient and C is the ratio of the 
velocity to the maximum velocity. The ratio, C, is assumed 
to be approximately 1 everywhere in the plume. The thrust 
coefficient Cp was evaluated using Eq. (4.33) of Shapiro,^ 
T is the isentropic exponent of the gas. 

For r constant, Hill and Draper^ give the following 
equation for the density as a function of 9 


-a 2 (i-cos e ) 2 

P = P0=O e 

where A is given by Eq. (3) . 


(4) 


It is assumed that the transition from continuum to 
free molecular flow occurs when the Knuds en number reaches 
one. For this study the Knuds en number is defined as 

g _ mean free path _ l_ ( 5 ) 

r r 

Using l = 3.736X10"25 (M/p&^) cm., (Ref. 5) where M is 
the molecular weight and 6 is the molecular diameter, we 
have by solving for p 

p = 3. 736X10" 25 -JL, _2_ ( 6 ) 

K 5 2 r cm 3 

Equating ( 6 ) to (4) and using (l) for p@ = o 

BK3 2 p r * 2 ,2 / ,2 

r = 1.0708X10 25 — £. e" A ( 1 - cos ®) (7) 

M 


which is the equation of the constant Knuds en surface in 
the plume . 


On the constant Knuds en surface the gas is assumed 
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Maxwellian with a mean velocity V in the radial direc- 
tion. From Holler's^ paper the density at a point 
P(x,y,z) due to molecules leaving an elemental volume dT 
adjacent to the constant Knudsen surface is 


fhL = *- 3 / 2 e - u2 


[ — u cos f + (i. 

2 V 2 


+ u 2 cos 2 e u2 cos2 f 



+ 


erf (u 


cos 


tfj 


m 


( 8 ) 


where f is the angle between the direction of V and the 
direction of flight of the molecules contributing to the 
density at P(x,y,z), u = | V| / (2RT/M ) , R the gas constant, 
T the temperature, and dQ is the differential solid 
angle of dT as seen from P(x,y, z). Using the result 
from another paper' 


d Q. = -2- sin 9 cos cp (l + A 2 cos 6 


- A 2 cos 2 d) - — + — (cos 9 - A 2 sin 2 9 


+ A 2 sin 2 9 cos 9) 


cos t 


sin 9 d 0 dcp 


® JL 


L \r 


(] (4 (f) 

- — + — cos e\ 

r e r e / 


sin 9 cos 9 


( 9 ) 

( 10 ) 


L 


2 =fc) 2 - 2 (i)(^) sl 


sin 6 cos 9 


>fjl\ (z-\ cos 0 + f-L \ 2 + U 
\ r e) \ T e) \ r eJ V r e 


J 


( 11 ) 


where cp is the angle that the projection of r in the 
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xy -plane makes with the x-axis . Note that r, 0, cp are 
the spherical coordinates. 

Equation (8) was integrated numerically over the top 
half of the constant Knuds en number surface. The upper 
limit for 0 was chosen as 0 p the maximum Prandtl- 

Meyer expansion angle from a Mach number of 5. The limits 
on <p were ±m/2. If r > r g when 8 = 9 p may - the con- 
stant Knuds en surface was assumed conical in shape with 
the cone half -angle being the expansion angle from the 
exit Mach number of 5 to the Mach number where the density 
is such that the Knuds en number equals one. For this 
region' 

cos 9 p cos cp - ^ sin e p )(f)(;r) sin e p a? 

( 12 ) 



where 0 is the Prandtl -Meyer expansion angle from Mach 
5 to the^Mach number such that K = 1 and L is the dis- 
tance from dx to point p. From Shapiro 


a p ■ tan ' 1 '\^TT (l l - - tan ' 1 


( C 


- |COS"“ si- - cos 


% 


where Mg, the exit Mach number, equals, 5. 
number where K = 1. 


Vrrr K - n 

(13) 

M e/ 

Mg is the Mach 


4 - rhffiT 1 ' ] 


For any given r < r e 
using K = 1. 


at 0 


p,max> 


(14) 

p is found from Eq. (6) 


The results for points in a plane located at the noz- 
zle exit are shown in Figs, 
and 
AA l 

a/ aJ- Ilu-?. fo 
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j-iits rcoui.ua xui jjuiaut) xij. a. j.ucaocu cx u uijlc iiu^j — 

exit are shown in Figs. 1, 2, and 3 for helium, argon, 
nitrogen, respectively. 
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Figure 3. 
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